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Abstract 

In this paper, we study several issues in the linear equation-of-motion (EoM) and in-in 
approaches of computing the two-point correlation functions in multi-held inhation. We 
prove the equivalence between this EoM approach and the hrst-principle in-in formalism. We 
check this equivalence using several explicit examples, including cases with scale-invariant 
corrections and scale-dependent features. Motivated by the explicit proof, we show that the 
usual procedures in these approaches can be extended and applied to some interesting model 
categories beyond what has been studied in the literature so far. These include the density 
perturbations with strong couplings and correlated multi-held initial states. 



1 Introduction 


The cosmic inflation 011 is a unique probe of high energy physics, providing access to 
energy scales that is inaccessible in colliders in the foreseeable future. The most important 
and well measured observable from inflation is the power spectrum of the primordial density 
fluctuations. The current observations can be well fit by a nearly scale invariant power 
spectrum with a slightly red tilt [^[^. The simple form of the primordial power spectrum 
will be tested or challenged by many on-going and future experiments. For example, future 
ground based telescopes can constrain the CMB polarization with a precision close to cosmic 
variance; future 21 cm experiments may push the measurement of primordial power spectrum 
to unprecedented precision [^, where the limitation of cosmic variance is as small as a relative 
error of 

The exciting future of cosmological experiments poses challenges for the computational 
methods of inflationary perturbations. The best established method to calculate the in¬ 
flationary perturbations is the in-in formalism [9 11 . It is an operator approach derived 


from first-principles, and is defined to compute the correlation functions to all orders in per¬ 
turbations. However, sometimes in practice, it is increasingly hard to compute high order 
corrections in the in-in formalism, because of the algebraic complexity, the spontaneously 
breaking of Lorentz invariance, and the subtleties of UV and IR divergences 12,13. More¬ 


over, non-perturbatively, the in-in formalism is not well defined, and in some situations 
provides formal expressions at most. 

Alternatively, solving the equations of motion (EoM) for fields is another way to compute 
the primordial perturbations. In contrast to the in-in formalism, the method of computing 
correlation functions to all orders in the EoM approach is not known, although there is a 
proposal on how to generate all the tree-level results [^. The method of linearly computing 
the two-point correlation functions in the multi-field inflation models, however, is well known 
in the literature. See e.g. 15 for earlier works, [Tbl for a review, and 17,35 for recent 


examples. Despite of the limitation, in the EoM approach, we deal with differential equations 
of c-numbers. This is particularly useful for numerical computations. In addition, as we will 
show, this method is even defined non-perturbatively at the linear level. 

The EoM approach used to linearly compute the two-point correlation functions has 
several special but well-known steps. Despite the general belief, an explicit and general proof 
for the equivalence between this EoM approach and the in-in formalism is not present. The 
equivalence has only been demonstrated in some special examples, e.g 17,19 . In this paper 


we provide such a proof. Inspired by the explicit proof, we also discuss several extensions 
beyond the known formalisms. 

This paper is organized as follows. In Sec. we first review the procedures used to 
compute the power spectrum in multi-field inflation models in both the in-in formalism and 
the EoM approach. In Sec. We prove the equivalence between the EoM approach and the 
in-in formalism at the linear order. In Sec. we use explicit multi-field inflationary models 
to check this equivalence. Motivated by the explicit proof of the equivalence in Sec. in 
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Sec. I^we show that the usual procedures in these approaches can be extended and applied 
to some new categories of model examples. These include the cases where the couplings are 
non-perturbative and where the initial states are correlated multiple helds. We conclude in 
Sec. m 


2 Tree-level power spectrum in multi-field inflation models 

In this section, we outline the two approaches of computing the two-point correlation func¬ 


tions in multi-held inhation models, namely linearly solving the EoM (Sec. 2.1) and com¬ 


puting tree-level diagrams in the in-in formalism (Sec. 2.2). We emphasize some seemingly 
ad hoc procedures required in the EoM approach. 


2.1 EoM approach 

Consider a multiheld inhation model with N number of helds, (pai'^Tt), a = 1,N. Denote 
(j)ai^,t) and as the background solution, and (50a(x,f) and S(pa{'^,t) as the perturba¬ 

tions. Denote L as the part of Lagrangian that is quadratic or higher in perturbations (50a 
and (50a, 


L{^8(j)ai ^0a) — -h(0a) 0a) -h(0a) 0a) 


'^3 dL 
a x-^d(Pa 
50a 


The EoM is 


d dL 


dt \d6^a 


dL 

86 (pa 


= 0, a = l,...,iV. 


81 

d^x^^6(pa . ( 2 . 1 ) 

50 a 


( 2 . 2 ) 


Here the partial derivatives are the functional derivatives. To linearly compute the two- 


point correlation functions, we only keep the quadratic terms in L. So (2.2|) is N number of 
coupled second order linear diherential equations. The mode function o: 
is dehned in the momentum space. 


50a, denoted as Ua, 


Ua{k,t) = / d^x5(pa{x,t)e 


2kx 


(2.3) 


These mode functions satisfy the same EoM (2.2) in the momentum space 


In this section, we consider the usual cases where the helds are decoupled initially. Con¬ 
sequently, each held takes the following initial conditions at the initial time to, 

Ma(k,fo) = Mr(k) , Ma(k,to) = <“(k) , o = 1, . . . , iV , (2.4) 

where the u™ and u™' satisfy the initial commutation condition. For example, if the helds 
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have canonical kinetic terms, these conditions are 


a^(to) {uaU*^ — c.c.) = i , (no snm over a) 


(2.5) 


We make a few clarification remarks on the terminologies “the initial commutation relation” 
used here and the “Wronskian condition”. Wronskian condition is a time dependent condi¬ 
tion for a decoupled single held mode, while the initial commutation relation is specihed only 
at the initial time to- If th® fields are decoupled initially, the initial commutation relation is 
the Wronskian condition evaluated at to- These conditions apply to both the Bunch-Davies 
vacuum and non-Bunch-Davies vacua, therefore in (2.4) the values of v}™ and u™' can be 


chosen arbitrarily as long as (2.5) holds. After to, we turn on interaction, so (2.5) is only 


needed at to but in general may not be satisfied for all t. 

The procedure of getting the two-point correlation functions is as follows. We need to 


solve the EoM (2.2) for N times separately with different sets of initial conditions. 


For example, for the 1st time, the initial condition is taken to be 


,(i) 


u\ = u 


.-.(b 


1 ) 


u\ = u 


1 ) 


u 


w = 0 , «!.’> =0 (a # 1) . 


(2.6) 


This leads to the 1st solution at the final time f, which we denote as u a\k,t), a = l,...,N. 
For the 2nd time, the initial condition is taken to be 


,( 2 ) 


Uo = U 




2 ) 


Ur, = U 


2 ) 


U 


<=> = 0 , «<=> =0 (a # 2) . 


(2.7) 


( 2 ) 

This leads to the 2nd solution at the final time t, Ua {k,t), a = 1,..., N. Repeat this for N 
times, and finally we get N sets of solutions at f, t), where a,a = 1,... ,N. 

The two-point correlation functions of the fields are determined by these solutions. 


N 


(50ap(t)(50bq(t)) = (27r)^^4“)(p,t)4“^*(q,f)^^(p 




( 2 . 8 ) 


a=l 


Notice that, if we directly supply all the initial conditions (2.4) to the FoM (2.2) and 


solve it only for once, we in general would have got the wrong results, even though the 
initial conditions are chosen such that the multiple fields are decoupled and the Lagrangian 
contains no interaction initially. 


2.2 In-in formalism 

Similarly, we denote 0a(x, t) and 7ra(x, t) as the background solutions, and 6(j)ai^,t) and 
67ia{'x,t) as the perturbations, where ^^(x, t) and 67iai^,t) are conjugate momenta associ¬ 
ated with the fields (j)a{'^,t) and S(j)ai^,t), respectively. We also denote H as the part of 
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Hamiltonian that is quadratic or higher in perturbations Scpa and 


H [(50„, 5Tla, t]= H [(pa, TTa] “ H [(pa, ^ j d^X^6(pa j d 


'^X^dTTa . (2.9) 

OTla 


To get the tree-level two-point correlation functions, we only keep the quadratic terms in H. 
The evolution of the fields is governed by the following EoM, 


, S'k^{x,t)^i H,STTa{x,t) 


( 2 . 10 ) 


a = 1,..., N. It is easy to show 10 that this EoM is equivalent to (2.2) 


The in-in formalism is a perturbative method of solving the EoM (2.10). We split if into 
two parts, 


// (5ir„, (] = ffo IS't'a, Sx.,, t] + H,lS(j>„,SXa,t] , 


( 2 . 11 ) 


where Hq is the quadratic kinematic part of the Hamiltonian describing the leading order N 
number non-interacting fields. The helds with this leading order behavior, denoted as S(pa 
and Snp, are called the interaction picture fields and follow from 




( 2 . 12 ) 


The initial conditions for this EoM is given by (2.4) and (2.5). Because interaction is not 


introduced in ffg, the mode function of each interaction picture field satisfies the Wronskian 
condition. 

The correlation function (Q) is then given by perturbatively expanding the following 
expression in terms of powers of Hj, 


(Q) = (^^1 


T exp ( i / Hi{t)dt 

/ to 


Q^it) 


T exp y—i [ Hi(t)dt 


'to 


in). 


(2.13) 


where 


Hi(t) = H,[S(pi,5xl,i\ , 
Q‘{t)=Q[S(,lSx‘,] . 


(2.14) 

(2.15) 


In the perturbative expansion, we plug in solutions of (2.12) in terms of the following de¬ 
composition. 


d(pa ^a®a,k T k i 

d'Xa = Waaa,-k + -k i ’ ( 2 - 16 ) 

where a = 1,..., iV labels different helds. Note that each held is expanded into one set of 
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the creation and annihilation operators^ 

The in-in formalism is derived from the hrst-principles [9 11, while the status of the EoM 
approach is not entirely clear given the procedures. For single held inhation, the equivalence 
of the above two approaches is obvious since they are solving the equivalent EoM, (2.2) and 
(2.10). For multiheld inhation, it is not obvious why we should solve the EoM separately for 
N times, instead of solving it once with all the initial conditions, even though in the model 
setup we are considering the cases where the initial helds are decoupled and non-interacting. 
In the following section, we justify this procedure in the EoM approach by showing that it 
satishes the same set of hrst-principles as in the in-in formalism. 


3 The proof of equivalence 


For completeness, let us start by briehy repeating the setup using a more detailed notation 
following |9 11 . Denote the Hamiltonian of the system as 


H[(j){t),7r{t)] = / d^x'H[0a(x,t),7rb(x,f)] , 


(3.1) 


where a,b = label the number of helds. We perturb (j)a and tTq around a time- 

dependent background, 4>a{'x,t) and 7fa(x, f), 

0a(x, t) = 0a(x, t) (50a(x, t) , 7ra(x, t) = 7fa(x, t) + STlai^, t) . (3.2) 


The background satishes the classical equations of motion, 

j , , dU ^ , dU 

0a(x,f) = — , 7ra(x,f) = -^. (3.3) 

OT^a OCpa 

We expand the Hamiltonian as 

/ fi-ij _ r fi-ij 

+ H[5(l){t),5n{t)]t] , (3.4) 

where we use H to denote terms of quadratic and higher orders in perturbations. 

If we impose the canonical commutation relation for <pa and tTq, we immediately get the 
same relation for 6(f)a and S^a, 

[50a(x, t), 67lb(j, f)] = iSabS^i^ - y) , 

[50a(x, t), S(j)b{y, t)] = [57ra(x, t),S7rb{y, t)] = 0 . (3.5) 


Using these relations, we can derive the equations that govern the evolution of the pertur- 

. In Sec. 


^This is the decomposition that has been used so far [9 
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5.2 


we discuss the extension. 
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bations 


<50a(x,t) 


H [6(j){t), 67r{t)-,t], 6(j)a{y^,t) 


57r„(x,t) 


H[6(j){t),67r{t)-,t],67ra{yi,t) . 

(3.6) 


However, for our purpose, let us not assume (3.5) for the moment. Instead we define 


N 

a=l 

a=l 


(3.7) 


where OQ^k,_k (o^ = 1) ■■•5 are N sets of the usual creation and annihilation operators, 
and satisfy the usual commutation relations. 


[a«,p, 4 ,_q] = ( 27 r)^ 5 ^(p + , 

K,p, a/3,-q] = 0 , 4 _ J = 0 . (3.8) 


Notice the difference between (3.7) and (2.16) is that in (3.7) each held is decomposed into 
N sets of creation and annihilation operators. 


As the initial conditions, we require that (3.7) satisfy the canonical commutation relation 


(3.5) at to, namely 


[50a(x,to),^vr;,(y,to)] =f5ab5^(x-y) , 

[(50a(x, to), 50fe(y, to)] = [5vra(x, to), 57rfe(y, to)] = 0 , (3.9) 


although we still need to show whether (3.5) are satished by (3.7) for all t. This initial 
condition translates into a number of conditions for the mode functions and w^\ 

To give an example without cluttered notation, let us look at the two-held model. The 
generalization is trivial. In the two-held model, we need two sets of creation and annihilation 
operators Ok, and &k, 


S4>1 = MiOk + Uifek + C.C. , 

S(j)2 = M2Cik + + C.C. , 

b'Kl = PlOk + Q'l&k + C.C. , 

57^2 = P2ak + <?2&k + C.C. . (3.10) 


As mentioned, the initial state usually studied are the decoupled initial state, which means 
that, at the initial time to, the mode functions in (3.10) are diagonal. 


ui = M2 = 0 , gi = p2 = 0 . 


(3.11) 
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In such cases, the mode functions Ui and pi, as well as V 2 and ^ 2 , satisfy the initial commu¬ 
tation condition, 


MiPi - u^pi = i , V2q2 - = I 


(3.12) 


Instead of deriving Eq. (3.6) from (3.3)-(3.5) as done in the in-in formalism [9 11 , here 


we start with this equations but define the commutators in Eq. (3.6) using the dehnitions 


in (3.8). At this moment it is not clear whether the commutation relations (3.5) and (3.8) 
are equivalent. Our next step is to construct solutions for v!a'^ and to achieve this 
equivalence. 

With this d ehni tion of commutation relations, we can determine the evolution of 
and using (3.6), which have the following solutions. 


6(j)a{^fi) = U , (57ra(x,f) = f/ {tfio)6Tra{^fio)U{t,to) , (3.13) 


where U satishes 

= -iH[6((){to),67i{to);t]U{t,to) (3.14) 

with the initial condition 


U{to,to) — 1 . 


(3.16) 


As mentioned, the initial condition 500(x, to) and 57ra(x, to) are chosen such that (3.9) is 


satished. Because of (3.13), once (3.9) is satished at to, the commutation relations (3.5) will 


be satished for all t. Therefore, we have ensured the equivalence between the commutation 


relations (3.5) and (3.8) for all t. The equations (3.6) dehned in two approaches are then 
equivalent. 


To summarize, what we have shown is that the decomposition (3.7) and the corresponding 


mode functions satisfy the same hrst-principle quantization conditions and equations of 
motion as in the in-in formalism, and therefore are equivalent to the in-in formalism in the 
validity regimes of the two methods. Now it is easy to see that this decomposition leads to 


the procedure in Sec. 2.1 Plugging (3.7) in (3.6) and matching the coefficients of the N sets 


of creation and annihilation operators, we can see that we effectively have the same EoM, 
satished by the mode functions m", for each set of creation and annihilation operators. The 
only difference is the initial conditions. For the decoupled initial state, the condition given 


by (3.11) and (3.12) reproduces exactly the condition given in the procedure in Sec. 2.1 


(2.5), (2.6) and (2.7), each time we solve the EoM. Note that to match both sides of (3.6) 


using (3.7), the EoM has to be linear. Using (3.7) it is also straightforward to see that the 


two-point correlation functions are indeed given by (2.8). 


The explicit proof also motivates the following considerations. 


In the in-in formalism, the perturbative expansion relies on the fact that the bilinear 
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interaction terms in Hj are perturbative with respective to Hq. On the other hand, 
since all we have shown is that the EoM procedure satishes the same hrst-principle, 
this procedure still applies even if Hj is non-perturbative, in which case the linearly 
coupled EoM can be solved at least numerically. 


The most general initial conditions that have to be satished are (3.9). The diagonalized 
initial commutation relations (3.11) and (3.12), which correspond to the case where 
the helds are initially decoupled, is only a special case. There can be more general 
cases where the helds are correlated initially. In the EoM approach, this corresponds 
to non-zero off-diagonal terms in the initial conditions for the mode functions in (3.7), 
and the rest of the procedure is the same. The in-in formalism can also be modihed 
accordingly. 


We study these two cases in more details in Sec. 


4 Examples 

In this section we explicitly check this equivalence using two non-trivial examples, namely 
the quasi-single held (QSF) inhation model and the standard clock (SC) model. In the hrst 
example, the potential has an approximate shift symmetry; in the second example, it contains 
some features. In each model, we study the cases with diherent vacuum choices, resulting in 
either scale-invariant corrections or scale-dependent features in density perturbations. 


4.1 Quasi-single field inflation model 


The QSF inflation model ^ 13,20^ 28 describes a two-held inflation model in which the 
inflaton is coupled to a massive held with mass of order H. The following is the matter 


sector Lagrangian for such an example 12,13 


C = -^{R + afg^’^d.ed^e - V,,{e) - - K(a) , (4.1) 

and we assume the following potentials for two helds 

Kr(0) = , (4.2) 

I4(cr) = ^mla^ . (4.3) 


See [^,13,20-28 for motivations and generalizations of such models. The held a has a mass 
comparable to the Hubble parameter and 6 is the light held driving inflation in a constant 
turning trajectory. The background evolution in the case of constant turning {& = 0) is 





















given by the following equations 


= ^{R + (Tof 9l + Kr (0o) + K (^^o) , 


+ {R + a,Y 
K (<^0) = (-R + CTo) 6l , 


= 0 , 


^ “ 7/2 


H {R + ao)H^ 


2MlH^ 


(4.4) 

(4.5) 

(4.6) 

(4.7) 


Perturbing the Lagrangian up to second order gives the following Lagrangian for the pertur¬ 
bations 




\{R + <^of 
2a^ReQ5a5e , 


59 --{d,59f 




(4.8) 


in which we have neglected the mass term for 59 and the second line is the bilinear coupling 
between two helds. Since cxo = const., we redehned R such that cxo = 0 in the above equation 
and hereaftei0 The full equations of motion for the two fields are given by 


59 + ?>H59 + —59 + 2 


3H9o + 9o 


R 


5a -|- 2,-^5a = 0 , 
R 


5a + 3H5a -|- 


P' 

V” + — 


5a - 2R9o59 = 0 . 


(4.9) 

(4.10) 


As we discussed in the previous sections our proof is not based on the BD initial condi¬ 
tions; so we can assume a more general, but still decoupled, initial condition as follows 


Ht 


Ruk = {Ce e-*"" - De e*"") , n, = 0 , 


and 




Ruk = 0 , Vk = i 


(4.11) 


Ht 


{Co 


—ikr 




(4.12) 


where Uk and Vk are the initial mode functions for 59 and 5a and the above initial conditions 
have to be set at sufficiently early times, i.e. /c|r| 3> 1. Initial conditions have to satisfy the 
commutation relations which are equivalent to requiring 


|C«P-|D»P = 1, = 


(4.13) 


^In the numerical computations we do not neglect the time dependence of cto as it may have some evolution 
at initial times. However, we do neglect the mass term of 59. 
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Using the above initial conditi ons o ne can solve the equat ions o f motion numerically twice 

The initial condition (4.11) leads to u^p\t), and (4.12) 


2.1 


following the procedure in Sec. 
leads to Up (t). The two-point correlation function is given by 

(sepSe^) = + l“h(*)l^) '5’(p + q) ■ 


(4.14) 


On the other hand, for the in-in approach, we use the so-called commutator form because 
it is more convenient in terms of avoiding some spurious IR divergence in the intermediate 
steps 12,13 . For a generic interaction Hamiltonian Hj and at the leading order in pertur¬ 
bations we have 


rti 


{66^) D - / dh dt2{[Hi{h),[Hj{u)je{ty]]) . 


'to 


'to 


Assuming the interaction Hamiltonian in the following form 

Hj{t) = C{t)6a6e 


(4.16) 


(4.16) 


yields 


{Sey = -4 Re 


dti / dt 2 C{ti)C{t 2 ) v{t 2 ) u{t 2 ) v*{ti) u*{t) {u{ti)u*{t) - u{t)u*{ti)) 




'to 


(4.17) 


in which u and v are the free held mode functions for 69 and 5cr, respectively. For the QSF 
model we have Hi = —2a^R9oSa69^ 

In the following hgures we plot the fractional correction to the power spectrum of curva¬ 
ture perturbation. The curvature perturbation is related to the light held 9 via 


C = -j«. (4.18) 

In Figjljwe plot the results for the weak coupling limit in which the two approaches agree 
pretty well. In Fig|^ we show the time dependence of the results from a specihc mode in 
the BD case showing that the equivalence is always valid between two initial and hnal time 
slices. 


4.2 Standard clock model 


In this subsection, we consider an example in which the potential contain features, namely 
the standard clock models. The Lagrangian of the example of the standard clock model 
in 29,30 is very similar to that of the QSF inhation model. There are two diherences in 


the model setup. First, the mass of the cr-field is much larger than H in the standard clock 


^Transforming from Lagrangian to Hamiltonian, a small mass term for 5a will be induced by the bilinear 
term. We will neglect this extra term. 
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Figure 1: The correction to the power spectrum for the QSF model in the weak coupling limit. 
The left (right) panel is the results for the BD (a non-BD) initial conditions. For the non-BD 
case we have set Dg = Dfj = 1/2 and hence Cg = = -y/S/d. For other parameters we have set 

R = 0.8Mp, mo — lAH and m ~ 0.06H which correspond to the coupling {0/H) ~ 0.004. 



Figure 2: The time dependence of the correction to power spectrum for the weak coupling case. 
The horizontal axis is number of efolds counted from the initial time. The parameters are similar 
to ones in Fig{^ The mode crosses the horizon around n ~ 6. 
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model 23,31 , so asymptotically the inflaton rolls in the trough of this cr-potential valley and 


the effect of the massive held is negligible. Second, an earlier phase that describes how the 
inhaton settles down to the bottom of this valley is added, during which the massive held 
plays an important role in density perturbations. The motivation for this model is given 
in 


32,33 and summarized in 30 


So the Lagrangian is still given by (4.1) but we modify the potential to 


V; = 14o [l - exp(-aVo-J)] + 


a 


(4.19) 


The inhaton initially starts somewhere at the plateau of this potential, and then rolls toward 
the bottom of the valley, oscillates and settles down to an ehectively single-held inhation 
model. The transition occurs when a held reaches to the critical value a ~ cx/. The 


background equations can be obtained by varying the Lagrangian (4.1) which results in 

N 9 Ao 1.0 


3M|i/2 = 1 

9q -\- 3H9(\ -\- 


H 

^ " ~1p ~ 


29q&q 


+ 


' 2"° 
K;(^o) 


(i?-|-(To) (i?-|-cro)^ 

K- (^o) - (-fi + O-q) 9g = 
(R + cTo)^ 9^ + 


2MlH^ 


- (9o) + Kt (o'o) , 

(4.20) 

0 , 

(4.21) 

0, 

(4.22) 


(4.23) 


For the perturbations we have the free held Lagrangian similar to the one in the QSF model. 


i.e. the hrst line in (4.8) (again we neglect the mass term of 69). There are several bilinear 
terms for this model but for the purpose of this paper we only pick up the following one for 
simplicity 


A/:, 


+ Rfao9o6a69 

rl 


The EoM are 


69 + 


3H + 


2an 


R + (Jq 

(3 -|- e)0odo + ^'od'o -|- 

6(t + 3H6a + 


69 + —69 


29o&^ 


H{R + ao) 


+ 




9o(To 

H 

9o&o 

H 


fo _ = 0 

H 

{R + ao) ^69 = 0 


(4.24) 


(4.25) 


(4.26) 


Finally, to solve the above equations we impose the following generic initial conditions 


Ht 


{R + ao(r)) Uk = {Cg - Dg = 0 , 


(4.27) 
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APj 


APj 



Figure 3: The correction to the power spectrum for the SC model. The left (right) panel is the 
results for the BD (a non-BD) initial conditions. For the non-BD case we have set Dq = D„ = 1/2 
and hence Ce = Ca = y^ 5/4. For other parameters we have set R = 2.7Mp, mo ~ 0.7H and 
m ~ O.IH, aj ~ 0.05Mp, Vq — 2.7 x 10“®Mp and the initial conditions for background fields are 
tuned such that the transition occurs at n = 6. 

and 


{R + ao{T))uk = 0 , 


Vk 


. Ht 

‘71 




(4.28) 


The procedure of solving the EoM, and the formula for the in-in formalism approach, are 
the same as in Sec. 14.11 

In Figj^we plot the fractional correction to the power spectrum obtained by both ap¬ 
proaches for the two cases of BD and non-BD initial conditions]^ In Figj^ we show the 
time-dependence of one of the modes in the BD case. We see that two approaches match 
very wellj^ 


5 Several extensions 


The procedures in the two approaches reviewed in Sec. 2T and Sec. 2^ are commonly used 
in the literature. In this section, motivated by the explicit proof of the equivalence between 
the two approaches presented in Sec. we discuss several extensions of these procedures 
that may be applied to new categories of models. 


^The examples used here are different from the best-ht examples in 


29 


Here we used the large held 


examples where the effect of the 2nd held is more important, although this cases does not ht the data. Our 
purpose is to show the equivalence between the two approaches. 

®We observe that even for a large correction the two approaches match pretty well. It seems even for 
this seemingly non-perturbative case the higher order corrections in the in-in formalism are still negligible. 
We currently do not understand why this is the case, but we expect generally this would not be true (see 


Sec. 5.1). 


13 





































Figure 4: The time dependence of the correction to power spectrum for the SC model. The 
horizontal axis is number of efolds counted from the initial time. The parameters are similar to 
ones in Figj^ The mode crosses the horizon at n ~ 6.4. 

5.1 Density perturbations with strong coupling 

As we mentioned at the end of Sec. the validity of the EoM approach does not rely on the 
perturbative condition. Therefore, solving the linearly coupled differential equations in the 
EoM approach produces the full tree-level non-perturbative two-point correlation function. 
Formally in terms of the in-in formalism, this would correspond to a non-perturbative re¬ 
summation of all the tree-level diagrams for the two-point correlation function. 

We emphasize that this is of course not the full non-perturbative result because the loop 
diagrams are not included; nonetheless it is an interesting subset. Another limitation is that 
this procedure only applies to the tree-level two-point correlation functions due to the lin¬ 
earity condition in the EoM approach assumed in the proof in Sec. Possible generalization 
would be very interesting. 

To study this in an example, we note that, in the QSF inflation model, the coupling 
between the inflaton and the massive held a is of order d/H. In the previous examples, 
this coupling is taken to be small so the correction to the leading power spectrum is small. 
For large 9/H^ the perturbative expansion in the in-in formalism breaks down, while the 
numerical computation in the EoM approach is essentially unchanged. In Figj^we plot the 
results for such a strong coupling case. We see large mismatch between the result from the 
EoM approach and that from the hrst order term in the in-in approach. This shows that 
the higher order corrections in the in-in formalism are non-negligible. 

It would be interesting to see if the EoM may be solved analytically, or if all the tree-level 
diagrams in the in-in formalism may be re-summed. We leave these for future investigation. 

5.2 Correlated initial states 

In this subsection we consider a more general initial condition in which the two helds are 
correlated. As we already discussed, we only need to satisfy initial commutation relations 
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Figure 5: The correction to the power spectrum for the QSF model in the strong coupling limit. 
The left (right) panel is the results for the BD (non-BD) initial conditions. For the non-BD case 
we have set Dg = D(j = 1/2 and hence Cg = Other parameters are tuned as follows: 

R = 0.05Mp, mo — lAH and m ~ 0.18H which correspond to the coupling (O/H) ~ 8.64. In both 
figures, the blue lines are the EoM results, and the red dots are the first order result from the in-in 
formalism. 


and the proof works even if the helds are correlated initially. We show here the consistency 
of the two approaches in this general case in an explicit example. 

There may be variety of reasons that the helds are initially correlated, for example, due 
to some interactions in a pre-existing phase. See 34 for an example. Our approach here 
is phenomenological. We parameterize as many such possibilities as we can and do not 
attempt to construct concrete models. We vary these parameters and study their effects in 
the density perturbations. We also use both the EoM and in-in approaches to compute the 
same problems and demonstrate their equivalence. 

We start from the following two correlated helds and their associate conjugate momenta 


S(f) = Uia -|- Vib + c.c. 

(5.1) 

TT^ = pia + qib + c.c. 

(5.2) 

6a = U 2 a + V 2 b + c.c. 

(5.3) 

TTcr = P 2 a + q 2 b + c.c. 

(5.4) 


These are 8 complex parameters or equivalently 16 real parameters. We need to satisfy the 
following commutation relations 


[50, TT^] = {ui,pi} + {vi,qi} = i 

(5.5) 

[6a,7i„] = {U 2 ,P 2 } + {V 2 ,q 2 } = i 

(6.6) 

[50, 5cr] = {mi, U 2 } + {^ 1 ,^ 2 } = 0 

(5.7) 

[50,7r„] = {ui,p2} + {^ 1 ,^ 2 } = 0 

(5.8) 

[6a, 71^] = {u 2 ,pi} + {v 2 ,qi} = 0 

(5.9) 

= {pi,P2} + {qi,q2} = 0 

(5.10) 
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in which we dehned 


{p,q} =pq* -p*q. (5.11) 

Note that the above six equations are constraints for the imaginary part only, as the real part 
of the combination in the l.h.s vanishes automatically. Hence, in the most general case, we 
have 6 constraint for 16 parameters and as a result 16 — 6 = 10 real independent parameters. 

To go further and have more control on the parameters we can restrict ourselves to the 
case in which the following parameters satisfy the usual commutation relations, same as 
those in the initially uncorrelated case. 


{tii.Pi} = {K 2 W 2 } = i- 


(5.12) 


Plugging these into (5.5) and (5.6) yields 


= {m2,P2} = 0. 


(5.13) 


This splitting of equations gives two more constraints so it appears that we have 8 free real 
parameters left. However, after this splitting, one can check that one of the equations would 
be redundant and can be reproduced from the others. So we actually have 9 independent 
parameters. We can further parameterize by noting that (5.13) implies that the combinations 
Viql as well as U 2 P 2 are both real. We thus can dehne 


qi = RiVi (5.14) 

P2 = R 2 U 2 (5.15) 


in which Ri are real parameters. Note that this last step is just a new way of parameterization 
and does not restrict the independency of parameters. Similarly we can rewrite (5.7) and 
(S by 


U1U2 + V1V2 = 24 (|Mip + |n 2 p) 


(5.16) 


in which A is real and the factor in the bracket in the r.h.s makes A dimensionless. We will 
see soon that the strength of the initial correlation is controlled by A. After some algebra 
we can show that vi,qi,U 2 and p 2 can be expressed as secondary parameters in terms of 
Ui,pi,V 2 ,q 2 ,Ri,R 2 and A, 


Vi = 


U 2 = 


A{q2 - R2V2 ){\ui? + \V2?) 
PlUl + qlV2 - Ri\ui\^ - R2 \v2\^ ’ 
A{pi - RiUi){\ui\^ + \V2\^) 

PlU\ + qlV2 - Ri\ui\^ - R2 \v2\^ ' 


(5.17) 

(5.18) 


As for the parameters in Eq. (5.12) we can simply write them as a combination of positive 
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and negative freqnency of vacnnm in Minkowski space to make snre that they satisfy their 
own commntation relation as well. That is 


Ht 




(r 

V 2 ~ [Co 


—ikr 




(5.19) 

(5.20) 


and their corresponding conjngate momenta. To satisfy commntation relations we reqnire 

- \D^\^ = 1, lar - \D.? = 1. (5.21) 

To snmmarize we have nine real free qnantities parameterized as follows 


Cflij D<j}, Co, Doj Ri, i?25 


(5.22) 


The hrst fonr parameters are complex nnmbers for the negative and positive freqnency 


vacnnm which have to satisfy (5.21). The remaining three parameters are real and are 


responsible for correlation. Tnrning off these parameters tnrns off the correlation as well. 
Note that A is mnltiplied to all correlation terms. That is, setting A = 0 tnrns off all the 
correlation terms. 

So far we have parameterized a snbclass of models with correlated initial states, in the 
following we check the consistency of two approaches using the above parameterization. 

There is an interesting difference between the correlated and un-correlated initial condi¬ 
tions in the in-in expansion. In the case of correlated initial condition we need only one inter¬ 
action Hamiltonian to obtain non-zero two-point function while in the case of un-correlated 
we need two. This is a direct consequence of the fact that in the former case two helds are 
correlated (even in the absence of interaction) and hence their contraction is non-vanishing. 
As a result, if the strength of the interaction Hamiltonian is typically a, the leading order 
correction to the zeroth order power spectrum is at the order of a in the correlated case and 
at the order of in the un-correlated case. 

If the helds are initially correlated, the leading order correction to the power spectrum 
is given by the following expression 


{59^) D i{n\ dti[H{ti),S9{tf]\n) 


(5.23) 


= —4i Im 


dtiC{ti) [u2{ti)ul{t) +n2(H)vt(^)] [ui{ti)ul{t) + vi{ti)vl{t)] 


in which we assumed Hi(t) = C(t)Sa69. 

The common phase of the two helds, 50 and da, is not observable. At least for the case 
considered here, the relative phase between the two helds is not allowed. This is implied by 


the relations (5.17 5.18) where we have to rescale 50 and 6a by the same phase. As a result 
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Figure 6: The zeroth order (top) as well as the correction to the power spectrum due to the 
coupling (bottom) in QSF model in the case of initially correlated case. The parameters are the 
following: R = 0.5Mp, niQ ~ lAH, m ~ 0.016H, = D„ = 0, i?i = 0, i22 = 0 and A = 10“^. 

of this observation we can omit two independent but irrelevant parameters as overall phases 
for both fields. Hence we effectively have seven real independent parameters. 

Another important consequence of initial correlation is that the free held results, i.e. 
when the held coupling is set to zero, is also modihed as we need to solve the decoupled 
equations twice with diherent initial conditions and the hnal result is the combination of 
both solutions. Hence in the following hgures we also present the free held results besides 
the correction to the free held due to the held coupling. Figs. iEI and show the results for 
the QSF model in three diherent set of correlation parameters. 

Although the strength of the initial correlation is controlled by the parameters A, the 
overall behavior of the hnal results is also quite sensitive to Ri and R 2 as well. For example, 
if |i?i| ^ 1 and |i? 2 | -C 1, the power spectrum grows rapidly as a function of scale (Figj^. 
If |i?i| > 1 and |i? 2 | ^ 1, there are some large spikes at largest scales in the power spectrum 
(Figj^. Finally, if either (|i?i| > 1 and |i? 2 | ^ 1) or (|i?i| ^ 1 and |i? 2 | ^ 1) then the power 
spectrum is oscillatory around an averaged value. In this case the amplitude and frequency 
of oscillations seem to be unchanged throughout scales (Figj^. The above observations show 
that the initial correlation is a rich phenomenon and deserves further investigation which is 
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Figure 7: The same plots as in FigJ^but with Ri = 1, i ?2 = 1. Other parameters are as follows: 
R = 0.5Mp, mo — 1.477, m ~ 0.01677, = D^j = 0, and A = 10“^. 
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Figure 8: The same plots as in Figsj^andj^but with i2i = 1, i ?2 = 0. Other parameters are fixed 
by the following values: R = 0.5Mp, mo — lAH, m ~ 0.016i?, Dfp = = 0 and A = 0.1. 
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beyond the scope of this paper. We also note that the parameterization we have chosen is 
not the most general one and we have one real parameter less than the most general case. 

Higher order correlation functions with such initial conditions can also be computed in 
the in-in formalism by extending the procedure reviewed in Sec. 2^ As we have done here, 
we replace the decomposition (2.16) with (3.7). The contractions between different helds are 
then modihed accordingly with the dehnition of (3.7). With these modihcations, Eq. (2.13) 
can be straightforwardly applied to perturbatively compute not only the power spectrum as 
we did here, but also any higher order correlation functions. 


6 Conclusion and discussions 


To conclude, in this work we have shown the equivalence between the EoM and in-in approach 
for cosmological linear perturbations. A few non-trivial examples are provided. 

Inspired by the proof, we extended the formalisms of both approaches beyond the types 
of models that are usually considered in the literature. This deserves further study. 


1. Non-perturbative linear perturbation. We used the transfer vertex in QSF inflation 
as an example. In the parameter space that we have studied, the correction from the 
transfer vertex in the EoM approach is smaller than that of the hrst order result in 
the in-in approach. It is interesting to carry out a more complete survey of parame¬ 
ter space. More importantly, since the more characteristic property of QSF inflation 
model, i.e. the intermediate shape non-Gaussianities, shows up in the level of non¬ 
linear perturbation, it is interesting to seek a prescription that treats the in-in and 
EoM approaches complementarily: Use the EoM approach to determine the transfer 
vertex, and then combine with in-in to calculate the three-point interaction. We hope 
to study it in a future work. 

2. Multi-held inhation with correlated initial states. We have parameterized the initial 
state and studied simple examples among all possibilities. It is interesting to further 
explore the effects in different parameter space, and to relate the current study to 
inhationary model buildings and hgure out scenarios in which the initial correlation 
show up. 


It is also interesting to seek for a general prescription of EoM approach at the nonlinear 
level for multi-held inhation. See e.g. 35 for linear order. 
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